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In this paper, we consider a simplified error-correcting problem: for a fixed encoding process, to 
find a cascade connected quantum ciiannel such that the worst fidelity between the input and the 
output becomes maximum. With the use of the one-to-one parametrization of quantum channels, 
a procedure finding a suboptimal error-correcting channel based on a semidefinite programming is 
proposed. The effectiveness of our method is verified by an example of the bit-fiip channel decoding. 
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I. INTRODUCTION 



Quantum error correcting [l|, |3|, |J{ is surely a neces- 
sary technique to protect quantum states against deco- 
herence and unexpected noise in quantum computations 
or communications Q ■ An error-correcting procedure 
is composed of encoding and decoding processes; the for- 
mer is usuahy done by embedding an input in a higher- 
dimensional Hilbert space, e.g., a single qubit is encoded 
as 



a\<d)+b\l) 
= a|0c) + 6|l. 



a|000) e (C^)^ 



(1) 



The decoding process denoted by a recovery channel TZ 
is in practice implemented by a combination of unitary 
operations and classical measurements, which is generally 
represented by 

^£(|0,)(0,|)^7^£(|0c)(0c|). 

The superoperator £ represents the occurrence of errors. 
The procedure of error correcting formulated in is to 
expand the input space appropriately and to design the 
recovery channel TZ such that the worst fidelity between 
the input and the output becomes maximum, i.e., 



mm{(j)c\TZ£{\<l)c) {(t)c\)\4>c 

\<f><:} 



(2) 



Especially, a necessary and sufficient condition for the 
perfect error-correcting, TZS {\(j)c) {(f>c\) — 10c) (0c |, was 
given in 0- When the condition is not fulfilled, how- 
ever, any analytic ways to design the recovery channel 
have been unknown. 



The min-max problem (O is still hard to obtain a 
global optimal solution for even by means of numeri- 
cal methods. Actually, has replaced the problem Q 
by the simple maximization problem of a special case 
of Schumacher's entanglement fidelity jQ|. However, the 
criterion will be inappropriate for practical purposes be- 
cause it is not based on the worst input. 

In spite of the difficulty, this paper proposes a numeri- 
cal method to solve the original min-max problem ^ for 
a fixed encoding process. The key idea is the relaxation of 
the problem to a convex optimization one with linear ma- 
trix inequality (LMI) constraints, i.e., a semidefinite pro- 
gramming (SDP). The methodologies of SDP's have been 
used in many fields including control theory and even 
quantum physics. For example, ^ successfully applied 
the method of SDP to the test distinguishing entangled 
from separable quantum states. We also find 9] made use 
of an SDP to obtain the optimal quantum channel which 
approximates certain desired qubit transformations. 

The concrete derivation of the SDP is the following. 
Any quantum channel is one-to-one correspondent to a 
positive semidefinite matrix with a linear equality con- 
straint [Tolllllll^ . which directly concludes the convex- 
ity of the set of quantum channels Hence Eq. Q 
can be rewritten as a convex optimization problem with 
respect to the recovery channel. Moreover, by relaxing 
the set of inputs and applying the S-procedure 0, 13 , the 
LMI constraints are derived. 

The suboptimal recovery channel, which is obtained 
from the SDP derived above, guarantees error correction 
for the worst-case input, unlike 'F]. Also, for the bit-flip 
channel, it will show almost the same performance as 
that of a special error-correcting code, the majority-rule 
code. 
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II. PROBLEM FORMULATION 

Let Ti. and K, be n-dimensional Hilbert spaces of an en- 
larged input quantum state and the corresponding out- 
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put, respectively. For example, H = (C^)®^ in Eq. iQJ. 
In this paper we identify K, with 7i as they have the same 
dimension. We denote CiTi.) the set of all matrices on 7Y. 
A quantum state p belongs to SiTL) — {p € C{H) \ p^ = 
p > 0, Tt p = 1}. An input-output relation of a quan- 
tum state is represented by a quantum channel; an input 
quantum state p € SlTi.) is transformed to 



k=l 

where Ek € 'C(7i) is a matrix with the size n x n and 
represents channel properties which are often identical 
to the occurrence of some errors. Note that 



E 

fe=i 



ElEk 



has to hold in order that the output p' satisfies Tr p' = 1 . 
Let us call the set {Ek} the error channel. Note also 
that an element of the error channel Ek and its unitary 



transformation E',. 



J UjkEj, where U 



jk) 



unitary matrix, have the same input-output relation as 
follows: 



where p' is given by Eq. Q . The difference between the 
input p — |(/)c)(0c| and the output p' is quantified by the 
fidelity / := {(l)c\p'\(f>c), which is bounded by < / < 1. 
It is known that the fidelity becomes 1 if and only if 
p' — |0c)(</'c|- The minimization with respect to the input 
|(/)c) means that the performance of the recovery channel 
is evaluated in the worst case. The best performance 
= 1 is attained if and only if the following condition 
holds. 

Theorem 1 Q]- There exists a set {Rk\ satisfying 
F({i?fc}) = 1, if and only if 



(ic|£'|£'fc|4) = OL^k^u 



(6) 



holds for all i,j,k,£. 

In Eq. ©, ajk is a constant to be determined by 
choices of Ej and Ek, and Su is Kronecker's delta. A 
remarkable feature of the theorem is that condition © 
is written by using the properties of only the error chan- 
nel {Ek} and the code space C. However, Theorem 1 
does not tell us any way to search an optimal or even a 
good recovery channel when the condition (|SJ is not sat- 
isfied. The main contribution of this paper is to present 
a procedure to overcome this critical drawback. 



J2 EkpEl J2 E'kPE'kl G S{n). 



(3) 



fc=i 



fe=i 



We now explain the problem to be investigated. In this 
paper, for the sake of simplicity we fix the encoding pro- 
cess, or equivalently we fix 7i, and concentrate on the 
decoding process only. An encoded input p = \(pc){4>c\ is 
restricted into the code space given by 

L 

C:= [\cbc)eH I \cbc)^Y.^k\kc)} cn, 
fc=i 

where {|fcc)} is an orthonormal system in Ti. In exam- 
ple O, \lc) = 1 000) and 1 2c) = |111) are chosen. The 
decoding process is represented by the recovery channel 
{Rk} satisfying X^fcli ^l^fc = where the number of 
the elements Rk, M, needs not to be equal to that of the 
error channel {Ek}. The recovery channel {Rk} is just 
connected to the error channel {Ek}, and it is designed 
so that the output given by 



M m 



(4) 



k=l 3 = 1 



is as close to the input p — \(j)c){4>c\ as possible. It is no- 
table that the optimal set of matrices {Rk} is not unique 
even if they exist, because {Rk} has unitary freedom as 
seen in Eq. Q. The error-correcting problem, which was 
originally considered in Q , is addressed by 



max F{{Rk}) 
F{{Rk}) 



min {(t>c\p'\(t)c 

|0c)6C 



(5) 



III. ONE-TO-ONE PARAMETRIZATION OF 
QUANTUM CHANNELS 

It is known that any quantum channel from S{'H) to 
S{Ti.) is in one-to-one correspondence with a positive 
semidefinite matrix acting on Ti*"^. We first restate this 
fac t by a slightly different way from the conventional one 
in [if. 

Fixing an orthonormal basis {\i) |j)}, we denote a 
vector in H®^ |$)) ^ E"j=i-^ul«> \ We some- 
times write \i)\j) instead of \i) (g) \ j) for simplicity. Let us 
introduce a specific vector in H'^^ given by 



^|fc)®|fc)' 



(7) 



fe=i 



where the asterisk denotes complex conjugation of each 
element Then |e)) is independent of the selection of 
orthonormal basis: i.e.. 



^ \ak) <8) \ak)* = ^ \bk) (8) |6fc)* 



(8) 



holds for any orthonormal basis {|afc)} and {|^fc)}. We 
should remark that the above fine property is not satisfied 
if |e)) is defined via the conventional way 12j as |e)) = 
X]fc=i 1^) ^ 1^)- The vector |e)) also has an important 
property expressed as 



{A(g)I)\e)) ^ {I (E)A^)\e)), yAeC{H), 



(9) 



where T denotes the matrix transpose. Actually, expand- 
ing a matrix A £ C{TC) by using an orthonormal basis 
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{\k)} as A = J the left hand-side of Eq. © 

becomes 



Similarly, the right-hand side of Eq. is calculated as 

which implies the equality 0. 

Here we define a positive semidefinite matrix Xi G 
£(7i®^) associated with a quantum channel p' = 
Efe XkpXl as 

Xi ■.= J2iXk®I)\e)){{e\{xl®I). 

k 

It turns out that the trace-preserving condition 
J2k-^k-^k = I corresponds to (Tr-^ ^ id)Xi = /, and 
the map p — *■ p' is written in terms of JCi as 



p' = (id®TrH)f(/0p^)Xi 



(10) 



Here, id denotes the identity operator on Ti.. Conversely, 
it is known that any quantum channel can be represented 
by using a positive semidefinite matrix as the above form 
(|1(J|I [ij. That is, Eq. IIUII defines a one-to-one corre- 
spondence between a quantum channel and a positive 
semidefinite matrix on 7i®^. Actually, Xi is invariant 
under the unitary transformation of the structuring ma- 
trices X'l^ — J2j^jkXj, where U — (ujk) is a unitary 
matrix, i.e., 

X^^Y.^Xk®I)\e)){{e\ixl®I) 

k 

= J2iX',®I)\e)){{e\iX',^®I). (11) 

k 

We next introduce another matrix expression for a 
quantum channel, X2 G £(7i®^). The introduction of 
X2 in addition to Xi is essential from the computational 
viewpoint, which will be explained in the last part of Sec- 
tion V. Let us define a vector associated with a quantum 
state p e S{H) as 



(p®/)|e))eH« 



(12) 



which is obviously in one-to-one correspondence with p. 
In particular, the vector representation of a pure state 
p = 10) (01 is given by 



(i0)(0i®/)Eifc>®ifc>* = i'/'>® 



(13) 



since |e)) is independent of the selection of Mul- 
tiplying an input-output relation of a quantum chan- 
nel, p' = J2k-XkpXl, by |e)) from the right, we have 



{p''E)I)\e)) =Y,ki^kpXl(g)I)\e)). On account of Eq. ©, 
this relation becomes 

ip' /)|e)) - E(^fe ® ^*k)(P ® 



From the definition (I12II . this equation is described by 

|P'»-^2|P», 

where \p')) = {p' (g) I)\c)) and 



X2 ■.= Y,Xk(gX^ 



k- 



The trace-preserving condition is rewritten by 
{{e\X2 = ((e| Y^Xk® XI = ((e| ^ I ® XjX^ = ((e|. 



Similar to Xi, we see that X2 is also invariant under 
the unitary transformation of the structuring matrices 
X'f. — UjkXj. The matrix X2 is associated with Xi 
through the rearrangement of the elements as 

m*x2mr^m\*x,\m*. 

This relation is independent of the selection of {\k)} be- 
cause of the property ((SJ • Since the rearrangement map 
is obviously homeomorphism, the set of Xi is equivalent 
to that of X2. Denoting this relation by Xi — <i>(X2), 
the above discussions are summarized as follows. 

Theorem 2. Any input-output relation of a quantum 
state is written by using Eq. ifT^ as \p')) = X\p)), where 
the matrix X is included in the set 

X = {Xe C{n®^) I ^X) > 0, ((c|X = ((e| }. 

The transformation ^{X) is defined with respect to an 
orthonormal basis {\k)} as 

{m*x\k)\er^{m*^xmir. 

The set X is convex, and its dimension is ~ = 
- N. 

The convexity of X is obvious. We also note that for 
El e X and E2 & X, we have E1E2 € X and £^2-^1 G 
X, which reflect that the fact the cascade connection of 
quantum channels is also a quantum channel. 



IV. RECOVERY CHANNEL DESIGN BASED 
ON SEMIDEFINITE PROGRAMMING 

Let us rewrite the problem ||SJ) by using the vec- 
tor representation of a quantum state H12|l . As seen 
in Eq. H13|l . the pure input p = \(j)c){4>c\ is repre- 
sented by \p)) = \(j)c)\(t>c)* ■ The output through the er- 
ror channel E and the recovery channel R is given by 
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\p')) = RE\p)) = RE\(t)c)\(t>c)* ■ Then, the inner product 
of \pf) and \p'J) is calculated as 

W)) = {ci^c\{<l^c\*RE\(l>,)\(t>cY 

= ^|(0,|i?ji?fc|(/)c)P = (0c|p'|0c), 
j-k 

which is just the fidelity between the input and output. 
In view of this, the problem ^ is rewritten as 



max FiR) 

FiR) := min RE\(j)c)\(l)c 

\4>a)&C 



(14) 



The optimal recovery channel Ro of the above problem is 
uniquely determined, while in the original problem jSJl, 
the optimal set of the matrices {Rk} cannot be deter- 
mined uniquely due to the unitary freedom Q. 

Before describing the procedure to find a suboptimal 
recovery channel, we derive a necessary condition for the 
perfect recovery channel to exist. The element 
of E^E is calculated as 

{z,\{z,\*E^E\j,)\j,y 

= {^,\{^,r[Y,{ElE,)<E>{ElE,r]\Jc)\Jcr 

^J^\{^c\eIe,\j,)\'. 



When the perfect recovery channel exists, the relation 
{ic\El,Ei\jc) = OiktSij has to hold for all i,j,k,£ from 
Theorem 1, which leads to 

m{ic\*E^E\jc)\jc)* = \ake6,j\^ := aS,,. (15) 

k,e 

Hence a perfect recovery channel never exists under the 
condition {ic\{ic\* E^ E\jc)\jc)* 7^ ctSij, which is easy to 
check. 

Now we shall give the procedure to find a suboptimal 
recovery channel when the condition H15|l does not hold. 
First note that the original problem H14|l is equivalent to 
the following minimization problem: 

min e s.t. ((/.e|((/)c|*fi-E|0c)|0c)* > 1 - e, V|(/.c) eC. 
We then relax the problem into 



mm e s.t. 



\RE\ 



> 1-e, V|0)) eC(»C*, (16) 



where C* is defined as a linear subspace spanned by 
{|fcc)*}- The meaning of the relaxation is the follow- 
ing: any element va C ® C* C H.®'^ is always given by 

l'?^)) = Sjife=i '^ifebc)|^c)*, where the coefficient <^jk e 
C has no restriction. Therefore, the original input 

|'/'c)|0c)* = Y^j,k=\ '^i-^fcb'c)|^c)* is obviously included in 
C ® C* , which indicates that the input is allowed to be 



a linear operator of the form (/) = j,^-^ 0jfc|jc)(^c| in 
addition to the pure state p = |^c)(^c| in the relaxed 
problem (|16|l . The relaxed condition is equivalent to the 
following inequality: 



-(i?£;) + -(fi£;)t + (e-l)/ 



> 0, 



V|0)) eC®C*. (17) 

Now we utilize a famous formula named the 5- 
procedure [tIIi^: 

ix\X\x)) > 0, V|x)) e 5 :- {|^)) ^ | {{x\S\x)) > 0} 
^ 3r > s.t. X -tS >0. 

Note that X and S need not be positive. This formula 
says that we may find at least one positive number t 
such that X — tS > is satisfied instead of checking 
({x\X\x)) > for all \x)) & S. 

We apply the S-procedure to obtain an equivalent rela- 
tion to the relaxed condition H17|l . Assume S is a negative 
semidefinite matrix acting on (C ^C*)'^. Then, all vec- 
tors |(/))) e C(g)C* satisfy (((/)|S'|0)) = 0, and the condition 
(|17|l is equivalently transformed to 

3t>0, s.t. ^{RE) + ^{REy + {e~l)I-TS>0. 

Consequently, our relaxed problem is to minimize e sub- 
ject to the following LMIs: 

^{RE) + ^{REy + ie-l)I-rS>0, (18) 



$(fi) > 0, r > 0, 
((e|H= ((el. 



(19) 
(20) 



with changing S e C{{C (g) C*)-^). This is a typical SDP, 
which enables us to find a suboptimal solution Ro achiev- 
ing the maximum fidelity F(Ro). 



V. APPLICATION TO THE BIT FLIP 
CHANNEL DECODING 

Let us consider the bit-flip channel with flipping prob- 
ability p: 



T{p) = paxpax 



wheiep+q — 1 and CTa; — |0)(1|-|-|1)(0|. When any encod- 
ing and decoding processes are not performed, the mini- 
mum fidelity between a pure input state and the output is 
q. In this section, we consider two kinds of encoded input: 
one is = a|lc> +b\2c) := a|00> +&|11), and the other 
is = a\lc) + b\2c) := a|000) + 5|111). They are per- 
turbed by the error channels T®^ and T®^, respectively. 
For the first case, we will examine the gap between the 
overguaranteed fidelity due to the relaxation and original 
minimum fidelity. For the second, we compare our nu- 
merical procedure with a famous error-correcting strat- 
egy, the majority-rule code. 



5 



For the error channel T®^, the interaction operators 
are given by 

E3 ^ y/pql <^ a^, Ei=qI®I. 



The corresponding matrix E = X^Li Ek^E^e C{n^'^) 
is expressed as 



E 



■ qEi 


VMEs 


VpqE2 


pEi - 


VmEs 


qEi 


pEi 


VmE2 


VmE2 


pEi 


qEi 


VpqE3 


. pEi 


y/PqE2 


^/pqE3 


qEi _ 



We first get 

{U{lc\*E^E\2,)\2,)* = 3pV +/ 0, 

which violates the condition p5|l : thus, there does not 
exist a perfect recovery channel. Now C ^C* is spanned 
by {|1,)|1,)M1,)|2,)*,'|2,)|1,)M2,)|2,)*}, and hence we 
can take S G ^{{C ^ C*)^) as the following negative 
semidefinite matrix: 



-/' 




-h 






-h 




-I' . 



where li denotes the 4x4 identity matrix and /' :— 
diag{0, 1,1,0}. 

Consider the case p = 9/10. By using MATLAB LMI 
Toolbox, the set of LMFs (jTHJl, (Cni, and ^ is solved for 
the minimum error e = 0.196 with a suboptimal solution 

Ro = 






aBi2 + PBi3 


aBi3 + l3Bi2 


Bli + Bii ' 





1B22 


1B23 








1B32 


1B33 





Bli + Bii 


aBi2 + (3Bi3 


aBi3 + (3Bi2 






(21) 

where Bij is a 4 x 4 matrix whose element is 1 and 
others zero. The parameters are given by a = 0.28, /3 = 
0.09 and 7 — 0.22, and O denotes the 4x4 zero matrix. 
We here round all entries of Ro off to two decimal places. 
Note that Eq. (|21|l is one of the suboptimal recovery 
channels derived from the relaxed problem H16() . That 
is, Eq. H21|l guarantees the minimum error e = 0.196 for 
extra inputs which are taken outside the code space. 

Let us compute F{Ro), the actual fidelity for the above 
suboptimal recovery channel Ro, in order to check the 
gap between the optimal one and the suboptimal one 
obtained by the proposed method. An arbitrary input in 
C is described by 



,) = a\lc) + b\2c) = 



, \a\' + \b\' = l. 



For the error channel without any recovery, the fidelity 
between the input and output is 



{cl)c\{M*E\cl)c)\cl)c 



2p' 



> 0.01, 



where 5ft denotes the real part. This implies that the min- 
imum fidelity is 0.01. Also, the minimum fidelity without 
any encoding and decoding processes is (7 = 0.1. Thus we 
can assert that our recovery channel (|21|) works quite well 
because it certifies the fidelity 1 — e = 1 — 0.196 = 0.804 at 
least. Now the input-output fidelity for the error channel 
accompanied with the recovery channel H21|l is calculated 
by 

{cl>,\{cl)o.\*RoE\cl),)\cl)or 

= 0.56m + P^ + 2{q^ + O.lSpq) \\a*b\^ + ^{a*'^b^) 



> 0.56m + P 



0.860, 



where the equality holds when a = 1 or 6 = 1. This ac- 
tual minimum value is slightly larger than the guaranteed 
worst value 1 — e = 0.804. This difference appears for the 
following reason: the input \4>)) € C (8) C* corresponding 
to the worst case is selected outside the code space C for 
the relaxed problem 116|l , while the actual worst input is 
shown to be {(pc) — |lc) or |0c) = |2c) in C. However, 
we may conclude that the recovery channel (|21|l is not so 
conservative because the difference is not too large. 

The next subject is to compare our numerical pro- 
cedure with the (improved) majority rule code, e.g., 
|001) |0), |011) ^ |1) for < p < 1/2 and 
joOl) ^ |1), |011) ^ |0) for 1/2 < p < 1. An input state 
= a\0) + b\l) is encoded into = a|000) + 5|111) 
and passes through the triple bit-fiip channel T^'^. When 
< p < 1/2, the output is decoded into 



p' ^p''{3~2p)c7,\^){(j)\(7. 

It turns out that the fidelity satisfies 



2(1 



2p)|0)((/.|. 



F ■ — 



■IP 



= q^il + 2p) + 2p2(3 - 2p) [lapl^p + ^{a*H'') 
>9'(l + 2p), 

where the equality is attained when a = 1 or 6 = 1. 
Similarly, we have -Flnaj > ^^(3 — 2p) when 1/2 < 
p < 1. On the other hand, we shall solve the LMFs 
(|18|l . (|19|l . and (|20|l with the negative semidefinite ma- 
trix S — diagjO, — 1, — 1, . . .}, where zeros appear in the 
1st, 8th, 57th, and 64th entries, and the others are all 
— 1. Let us consider the case p = 9/10 again. Then the 
minimum fidelity via the majority-rule code is given by 
minFmaj = 0.972, whereas our numerical method yields 
a suboptimal recovery channel with the minimum error 
e — 0.048, or equivalently, the worst fidelity \ — e — 0.952, 
i.e.. 



\RoEV 



> 0.952, V|0)) GC(g)C* 
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The guaranteed fidelity is slightly less than min_Fmaj due 
to the relaxation. For the other case p = 1/10, anal- 
ogously, we have minFmaj — 0.972 while the LMI's are 
solved for the same worst fidelity 1 — e = 0.952. The 
above investigations conclude that the proposed method, 
which is applicable for general quantum channel without 
any prior knowledge, has almost the same performance 
as that of a special code for the bit-flip channel. 

We lastly remark on an important reason why R is in- 
troduced in addition to <&(-R). Actually, for the problem 
in 0, i.e., approximating certain desired qubit transfor- 
mations via a quantum channel, the solution does not 
depend on whether we use R to describe the LMI's or 
not. Now the fidelity to be optimized can be described 
in terms of $(-R) as 

F= min (0,1(0,1* ^(/®i?T)ci>(fi)(/®i?*)|^,)|0,)*, 

\<f>a)£C ^ 

and thus the same procedure as in Section IV leads to an 
LMI: 

^(/ ® El)<^{R){I ® El) + (e - 1)7 - t5 > 0. (22) 
fe 

However, the solution R!^ obtained from Eqs. p9(l and 
H22|l and (Tr-^ ® id)$(i?) = / is no longer the same as 
the one via the LMI's (HHJ), (O, and (HO)), because the 
relaxed constraints disagree each other, i.e., 

mY^ii ® Ei)'^{R){i ® Eim) ^ {{mm)- 

k 



Indeed, for the double bit flip channel T®^ with p = 9/10, 
the LMI's including Eq. H22(l provide a solution for the 
minimum error e = 0.749, whereas e = 0.196 via the 
LMI's lO, (Uni, and For this reason, unlike 0, 

the introduction of R in addition to $(-R) is essential in 
our problem formulation. 



VI. CONCLUSION 

In this paper, we have considered a simplifled error- 
correcting problem: for a fixed encoding process, to find a 
cascade-connected quantum channel such that the worst 
fidelity between the input and output becomes maxi- 
mum. With the use of the one-to-one parametrization 
of quantum channel 0, 0, ^2 > ^ suboptimal recovery 
channel can be determined as a solution of a semidef- 
inite programming. The effectiveness of the proposed 
method has been verified by studying the bit-flip chan- 
nel. Although we could find a suboptimal recovery chan- 
nel which is very close to the optimal one for the example, 
the condition where a good recovery channel exists has 
not been cleared yet, which is an important future work. 
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